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ABSTRACT 

We present a formalism describing the physical content of cross-correlation functions 
between a diffuse background and a population of discrete sources. The formalism is 
used to interpret cross-correlation signals between the unresolved X-ray background and 
a galaxy population resolved to high redshift in another spectral band. Specifically, we 
apply it to the so-called faint blue galaxy population and constrain their X-ray emissiv- 
ity and clustering properties. A model is presented which satisfies the recently measured 
constraints on all 3 correlation functions (galaxy /galaxy, background/background and 
galaxy /background). This model predicts that faint galaxies in the magnitude range 
B = 18 — 23 (covering redshifts z < 0.5) make up ~ 22% of the X-ray background 
in the 0.5 — 2 keV band. At the mean redshift of the galaxy sample, z = 0.26, the 
comoving volume emissivity is px ~ 6 — 9 x 10^® /i ergs s~"'^Mpc~^ . When extrapo- 
lated to fainter magnitudes, the faint blue galaxy population can account for most of 
the residual background at soft energy. We show how the measurement of the angular 
and zero-lag cross-correlation functions between increasingly faint galaxies and the X- 
ray background can in principle allow us to map the X-ray emissivity as a function of 
redshift. 
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NTRODUCTION 

o 

'^tiie X-ray background (XRB) has been the first diffuse cos- 
iMS; radiation ever detected over three decades ago (Giac- 
ni et al. 1962), but unlike the more famous microwave 
background discovered shortly after, its puzzling origin is 
yet to be solved. Although it is now known (e.g. Shanks et 
al. 1991) that quasars contribute a significant (> 30 %) frac- 
tion of the soft X-ray part of the radiation (~ 1 keV), the 
observed steep spectra of these objects, or of any others, do 
not resemble that of the hard X-ray emission. 

In recent years the XRB has been studied by means of 
analyzing the total intensity, the spectrum and the spatial 
fluctuations (e.g. Boldt 1987 and Fabian & Barcons 1992 
for reviews). In particular, the spatial fluctuations were an- 
alyzed by: (i) Source identification of high-flux regions (e.g. 
Shanks et al. 1991) (ii) Auto-correlation function (e.g. de 
Zotti et al. 1990, Carrera et al. 1993) (iii) Cross-correlation 
of the XRB with nearby galaxies (Turner & Geller 1980, 
Jahoda et al. 1991, Lahav et al. 1993, Miyaji et al. 1994, 
Carrera et al. 1995, Barcons et al. 1995). 

In this paper we extend the cross-correlation approach 
to higher redshift and confront it with recent observations 
(Roche et al. 1995). We present a general formalism which 
takes into account the clustering of the sources, their evo- 



lution, and the world geometry. We then apply this formal- 
ism to the so-called faint blue galaxy population. The un- 
expectedly numerous sources observed at faint blue magni- 
tude (Tyson 1988, Lilly et al. 1991, Metcalfe et al. 1992, 
Lilly 1993) have been a subject of active research in the 
last few years. It now seems that most of the 'excess' galax- 
ies (with respect to the volume density expected from local 
surveys), are narrow emission-line objects at moderate red- 
shift (^ < 1), with star-forming or Seyfert activity (Broad- 
hurst et al. 1992, Tresse et al. 1994). The occurrence of close 
pairs and disturbed morphologies among them (Glazebrook 
et al. 1995b, Driver et al. 1995) suggests that these object 
could be in the process of merging, thus explaining their 
much lower number density in local surveys. Or else, they 
may have declined in luminosity rather than number due 
the short time scale of their activity (e.g. Babul & Rees 
1992). Whatever the actual process of their apparent van- 
ishing, their starburst and/or Seyfert activity makes them 
likely X-ray sources, and therefore an XRB candidate popu- 
lation given their extremely large surface density on the sky 
(about 10^ objects per square degree at B ~ 23). 

In the next section we review the necessary ingredi- 
ents of cross-correlation analysis, namely the basic statisti- 
cal properties of the galaxies, of the XRB and of the source 
clustering. In the third section, we write the cross-correlation 
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and auto-correlation functions as a function of these ingredi- 
ents, both in Fourier space in terms of power-spectrum, and 
in real space using more specific assumptions and approx- 
imations. These functions have been explicitely derived in 
Appendix A. We then constrain the parameters (Section 4) 
using a recent measurement of the cross-correlation function 
between a magnitude limited sample of faint blue galaxies 
and the soft XRB (Roche et al. 1995), and present a model 
for the faint blue galaxy contribution to the XRB. We show 
that X-ray emitters (such as starburst or Seyfert galaxies) 
among the faint galaxies can account for the observed statis- 
tical properties of the XRB and propose further statistics to 
quantify their X-ray contribution as a function of redshift. 
Our conclusions and future work are discussed in section 5. 



2 BASIC PROPERTIES AND ASSUMPTIONS 

In this study, we cross-correlate known galaxies with the dif- 
fuse XRB. The calculations involve properties of the galax- 
ies, of the unknown X-ray sources, and their spatial cluster- 
ing. 



2.1 The galaxies 

Let us assume a population characterized by a luminosity 
function (f>i,(Li,, z) at frequency and redshift z, and a spec- 
tral energy distribution of the form oc (where /3 is 
the spectral index in the relevant band Ai^). In a survey 
with flux limits fmtn and fmax at frequency i^, the number 
density of resolved galaxies at redshift z will be: 



n{z) 



(1) 



where Lmin,maj:{z) = 47rr, {z){l + zy fmtn,max and ri{z) 
is the luminosity distance. We call Pi,{z) the total volume 
emissivity of the population observed in the band v ± Ai^, 
i.e. originating from the frequency band z)i'±{l-\- z)IS.i': 



pAz) = (1 + 2)' 



-13+1 



(pviLi,, z)Li,ALi,. 



(2) 



In the case of optically selected galaxies, the observables are: 
(i) the number of galaxies N{m, z)AmAz per deg^ on the sky 
at magnitude m and redshift z, measured from deep spec- 
troscopic surveys down to an optical magnitude of B Si 24 
(e.g. Glazebrook et al. 1994a), and (ii) the galaxy number 
counts C(m) = j N{m,z)Az, observed to B Si 27 (e.g. 
Tyson 1988). We use the following fit to the galaxy redshift 
distribution in the Bj-band proposed by Efstathiou (1995): 



N(m,z) = C(m) 



3z' 



2zi(m) 



exp 



Zc{m] 



3/2 



(3) 



where: 



Zc{m) 



0.0113(m - 17)^^^ -I- 0.0325 if 17 < m < 22 
0.0010(m - 17f -I- 0.0325 if m > 22. 



The evolution of the galaxy luminosity function results in 
a steepening of the faint end slope, i.e. an increase in the 
number or/and luminosity of locally faint and blue galaxies. 
The volume density n{z) of catalogued galaxies at redshift 
z (Eq. 1) appearing in the correlation functions verifies: 



Table 1 . Parameters on the X-ray background and X-ray sources: 





Soft band 

0.5-2 keV (ROSAT) 


Hard band 

2-10 keV (HEAOl) 


/ 


~ 1.3 - 2.1 X 10-** 


5.2 X 10"** ergs s~^cm~^sr~^ 


PxRB 


0.4 - 1.1 


0.4 


PO 


7 


~ 10^^ h ergs s-^Mpc-^ 


lAGN 


~ 3 


7 


Pagn 


~ 0.7 - 1 


~ 0.7 - 1 



n(z)Av(z) = Az I N{m,z)dm, 



(4) 



where Av(z) is the volume element per unit solid angle at 
redshift z (Appendix B). The mean projected number of 
galaxies per unit solid angle is: 



N = I n{z)<lv{z) : 



C(m)dm. 



(5) 



2.2 The X-ray background 



Resolved X-ray sources, essentially AGN's with ~ 1 keV- 
fluxes > 10~^^ ergs s~^cm~^, contribute ~ 30 % of the 
total extragalactic X-ray intensity (Hasinger 1992, Shanks 
et al. 1991). In this paper, the 'X-ray background' refers 
to the unresolved fraction of this radiation. Let us call 
Px{z) = Pd{z) + Ps{z) the volume emissivity of the XRB 
due to diffuse radiation and unresolved discrete sources at 
redshift z, observed in a given X-ray band. The diffuse com- 
ponent pd, first thought to be dominant and due to hot 
inter-galactic gas, has since been shown to be strongly lim- 
ited. Too much hot gas in the Universe would indeed distort 
the microwave background spectrum, and the COBE ex- 
periment has proved that no such effect was observed to a 
very high accuracy (Mather et al. 1990). The discrete source 
component ps is therefore strongly dominant, although no 
known X-ray source population has yet been shown fully 
consistent with the background properties, ps may be writen 
as the added contributions of various populations emitting 
in this band, i.e. a sum of Eq. 2's. In the following, we will 
assume the diffuse component of the XRB to be negligible 
{Pj;{z) Si ps{z)), and model the observed yi-ray light density 
of the sources in a given X-ray band as: 

Psiz) = poil + zy. (6) 
The mean intensity of the XRB per unit solid angle is: 

Px{z 



I 



4Trrf{z) 



dv{z) 



Po_ 
4n 



(1 + z 



,9-2 



Fdrc 



(7) 



where Fdrc = dv/r^ is a geometry dependent function of 
z defined in Appendix B. The soft (E ~ 1 — 2 keV) and 
hard (E > 2 keV) XRB exhibit vary distinct properties (see 
Boldt 1987 for a review) and we actually possess very dif- 
ferent data depending on the energy band. The known (but 
somewhat uncertain) parameters of the XRB and sources 
(based on Boldt 1987, Eabian & Barcons 1992, Hasinger et 
al. 1993, Boyle 1994, Chen et al. 1994) are summarized in 
Table 1. The measured intensity in the soft band actually 
varies within a factor of almost 2, depending (among other 



Me 



valu 



ill be indicated as either X or < X >. 
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things) on the Galactic component and resolved source re- 
moval procedures. In the following sections, we shall use 
T = i X 10~^ ergs s~^cm~^sr~^, where 1.3 ^ J ^ 2.1 in the 
0.5 — 2 keV energy band. 



2.3 Galaxies and X-ray sources 

We shall characterize the cross-correlation function between 
the hypothetical X-ray sources and the known galaxies by 
the usual function {(r, z). The excess probability of simulta- 
neously finding a galaxy and an X-ray source in the volume 
elements SVi and SV2 respectively, with separation ri2, is 
defined as: 



< n(SVi)ps(SV2) >= [1 + i(ri2,z)]n(z)SVips(z)SV2. 



(8) 



The spatial auto-correlation functions of nearby populations 
(e.g. galaxies and clusters) are well approximated by power- 
laws. For the cross-correlation function between the X-ray 
sources and the faint galaxies we shall assume a similar form, 
taking into account clustering evolution with redshift (Pee- 
bles 1980): 



{(r,^) = (l + ^)-(^+^' - 



(9) 



where r is the physical separation between the sources and e 
an arbitrary evolution parameter. Theoretical models com- 
monly assume 7 — 3 < e < for the galaxies. The upper 
limit corresponds to 'stable clustering' while in the lower 
limit case, the growth of density perturbations is due to 
the expansion of the Universe only. For local bright galax- 
ies, the correlation function parameters are 7 ~ 1.8 and 
ro ~ 5 ft~^Mpc (Peebles 1980). Faint galaxies on the other 
hand seem to be significantly less strongly clustered, with 
ro ~ 2 h-^Mpc (Ffstathiou 1995). 

To allow for non power-law models and more direct com- 
parison with theoretical models, it is convenient to express 
the cross-correlation in terms of power spectrum in Fourier 
space: 



P{k,z) 



m 

(l + z)' 



(10) 



where e' is an evolution parameter. (In the case where the 
power spectrum is a power-law with slope n, e' = e — n). 
As a working hypothesis, we assume all three correlation 
functions to be the same (as would be the case for instance 
if the faint galaxies were the only clustered X-ray sources), 
but the formalism can easily be modified to allow the auto 
and cross-correlation functions to be different. 
Operationally, the statistical angular galaxy /background 
cross-correlation is performed by dividing the sky into cells 
of small solid angle and calculating : 



r]gb{0) =< (N - N)e-{I - I)e-+e >ceHs, 
or, in a normalized form: 



Wgb{0) 



NI 



(11) 



(12) 



where the average is over all pairs of cells with angular sep- 
aration 9. 

We use the standard Robertson- Walker metric relations for 
the world geometry (distances, volume, etc.). These rela- 
tions are recalled in Appendix B. Unless otherwise stated. 



we assume an Q = 1, A = universe, and write the Hubble 
constant as Ho = 100 h km/s/Mpc. 



3 CORRELATION FUNCTIONS 

We have derived the formalism required to interpret cross- 
correlation and auto-correlation analysis. (Details of these 
somewhat cumbersome calculations are given in Appendix 
A to avoid overloading the main text). We first write the gen- 
eral formulae in Fourier space and then propose a more prac- 
tical version in real space, applied to the faint blue galaxy 
population. 



3.1 Angular correlation functions 

Given two functions of galaxy positions, X and Y, the angu- 
lar cross-correlation function between them at angular sep- 
aration 9 can be written as: 



r]XY{9) 



P(k,z) k gxY(k9) dk, 



(13) 



where P{k,z) is the power spectrum at redshift z and gxY 
is an appropriate window function of X and Y. For the 
galaxy /galaxy and background/background angular auto- 
correlation functions respectively: 



ggg{k9) = -!- / n^{z).Mkrc9)rtFdrc 

ZTT 



g,,{k9) = ^ j ] JoikrjytFdr, 



(14) 
(15) 



For the galaxy /background angular cross-correlation func- 
tion: 



9gb{k9) 



1 

2^ 



n{z)!^^Jo{krc9)rtFdrc 



(16) 



In the above and following equations, r;, fc, dv and F are 
functions of z defined in Appendix B. We omit to write the 
z dependency in order to make the equations more compact. 
Jo is a Bessel function of order zero, the definition of which 
appears in Appendix A. To illustrate the scales probed by 
the three statistical measures as a function of angular scale, 
we have plotted the three window functions gxY (assuming 
an arbitrary set of parameters described in the caption) for 
9 = 1' (upper three lines) and ^ = 1° (lower three lines). The 
thick solid and dashed lines represent P{k)k^ for a standard 
CDM model and a low density CDM model respectively. 
The product P(k)k^ x gxY is the contribution per logarith- 
mic interval d In A; centered on wave number k to the inte- 
gral in Fq. 13. On the arcminute scale, the kernels are con- 
stant with k for all three statistics down to very small scales 
{k~^ ~ 0.1 ft~^Mpc). As they are functions of k9, increasing 
9 simply results in shifting the kernels towards smaller k. For 
example on the degree scale, the three statistics are less sen- 
sitive to scales smaller than k~^ ~ 10 ft~^Mpc. Note that the 
XRB statistics is flux- weighted (oc r^"^) , hence less sensitive 
to larger scales compared with the number-weighted galaxy 
statistics. However, as in our model the XRB sources are 
distributed over a much larger range of redshifts, the back- 
ground auto-correlation probes larger wavelenghts than the 
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galaxy auto-correlation (and the galaxy /background cross- 
correlation probes intermediate scales). 

FIGURE 1 



3.2 Zero-lag cross-correlation function 

At zero-lag (9 = 0), we consider cells of solid angle to. A 
Poisson term, due to the potential X-ray emissivity pg of 
the catalogued galaxies themselves, adds to the clustering 
term: 



no = np + ic- 

The poisson term is defined as: 



rip = LoAIg = Lo 



4nrf 



dv. 



(17) 



(18) 



where A/g is the X-ray intensity contributed by the cata- 
logued galaxies. The clustering term can be writen: 



P{k,z) k go{k) dk, 



(19) 



where the kernel go is given by: 



27r 



3.3 Modeling 

With the various prescriptions and approximations de- 
scribed in Section 2 and Appendix A, the angular corre- 
lation functions result in simple power-laws of the angular 
separation 9 with power index 1 — 7. We can write all three 
functions in the form: 



nxY{0) = AxyO 

where: 



1-7 



(20) 



-4gg = H.yi / dz{l + Z 



A—1 



Pd.. \ I N(m,z)dm ) (21) 

A,, = H.,r2 [^^y I dzil + z f+'^^-'Ki-'F^ (22) 



ndfc 



A-gb — FlfTQ — — X 



47r 



dz{l + zf^'^-'r',-'' / N{m,z)dm (23) 

where p = 7 — (3-|-e) is a parameter describing the clustering 
properties of the galaxies with the XRB sources. Assuming 
for simplicity square cells with size uj = , the cluster- 
ing term (Eq. 19) in the expression of the zero-lag cross- 
correlation can be written (Appendix A): 



(24) 



we are interested in essentially are: the volume emissivity of 
the X-ray sources as a function of redshift, parameterized by 
po and g, and the clustering parameters ro and p. The auto- 
correlation functions of both the galaxies and the XRB have 



been studied elsewhere by various authors (e.g. Peebles 1980, 
de Zotti et al. 1990, Carrera et al. 1993). Our purpose here is 
to emphasize the cross-correlation techniques between high 
redshift galaxies and the diffuse background. The faint blue 
galaxies selected in a given magnitude range, stand as bea- 
cons allowing us to probe the XRB in the particular redshift 
range these galaxies are predominantly populating. Thus, 
measuring the amplitude Agb of the cross-correlation func- 
tion in various magnitude ranges can provide constraints on 
the clustering parameters and the volume emissivity of the 
X-ray sources. 



4 THE FAINT GALAXY CONTRIBUTION TO 
THE XRB 

4.1 Constraining the parameter space 

We have fitted the measured angular cross-correlation be- 
tween the XRB in the 0.5 — 2 keV band and faint galaxies 
in the magnitude range B = 18 — 23 (Roche et al. 1995) by 
the following power index and amplitude: 



(25) 



7 = 1.85 ±0.05 
A, = (1.5 ± 0.4) X 10"^ X IN 

The observed values for the mean intensity of the unresolved 
XRB in the above energy band and for the mean number 



of 18 < B < 23 galaxies are taken to be: / 



I X 10" 



ergs s~^cm~^sr~^ where 1.3 ^ i ^ 2.1, and = 2.6 x 
10^ galaxies per steradian respectively. Eor the purpose of 
exploring the parameter space, we introduce K = p + q — 2 = 
7 — e-|-g — 5asa global evolution parameter, and define the 
function: 

/. /.B = 23 

fm= (l + zfrl--' N{m,z)dm. (26) 

Jz J B=l& 

(We assume B = Bj -\- 0.2). With these new parameters, 
we can rewrite Eq. 23 in the following compact parametric 
form: 



Pori 



AttA^ 



(27) 



which we use to constrain the background emissivity and the 
clustering length. We consider a range of theoretical evolu- 
tion models for the clustering of the galaxies with the X-ray 
sources: —1.2 < e < 0, and a range of evolution for the light 
density of the X-ray sources: < g < 4 (from no-evolution 
to quasar-like evolution). These ranges imply —3 < K < 3, 
which we assume to be a reasonable and wide enough range 
for this parameter. Eigure 2a shows po as a function of K 
(Eq. 27). 

EIGURE 2a 

We find: 6 x 10^** < por^ < 17 x 10^** (where ro is in 
ft~^Mpc and po is in ih ergs s~^Mpc~^). Assuming 2 < 
ro < 6 ft~^Mpc and 1.8 < 7 < 1.9, this result implies: 



2.3 X 10" <po< 4.9 X 10^** ih ergs s"^Mpc" 



(28) 



It is clear that measuring the angular cross-correlation func- 
tion alone is insufficient to deduce the X-ray emission of the 
galaxies, unless very specific assumptions are made. Assum- 
ing as we did that the faint galaxies essentially cluster with 
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themselves, i.e. that the cross-correlation between the back- 
ground and the galaxies arises from their auto-correlation 
as opposed to their potential clustering with other X-ray 
sources, po parameterizes the galaxy own X-ray emissivity. 
In this case the background auto-correlation function (ACF) 
derived in Eq. 22 is the galaxy contribution to the XRB 
fluctuations. (Other clustered populations may produce ad- 
ditional fluctuations). Recent determinations (Chen et al. 

f994) ^ of the soft XRB anisotropy set the following 2-a 
upper limit: 



wtt{0) < 3.86 X fO" 



(29) 



on scales ^ ~ 2 — 300 arcminutes, or, in the present notation 
(Eqs. 20 and 22): 



Atb < 9.6 X fO'^r 



(30) 



Using both auto and cross-correlation measurements and 
combining Eqs. 22 and 27, we derive another constraint on 
Po (independent of ro ) as a function of the evolution param- 
eters: 



Po < 1.2 X fO" 



P /(p + g-2) 
J^{1 + z)P+^ii-^)rl-'' Fdrc 



(31) 



where / is the function defined in Eq. 26. The two solid lines 
in Eig. 2b bracket the above upper limit as a function of q 
when the parameter p = j — e — 3 varies from —1.2 (upper 
line) to (lower line). (The error bars on At are ignored in 
this calculation). 

EIGURE 2b 



4.2 The faint blue galaxy contribution 

The galaxy auto-correlation function rjgg has been measured 
by several groups. At faint magnitude, the data are consis- 
tent with low ro in contrast with the observation of bright 
local galaxies (Efstathiou 1995). We assume the fitting val- 
ues of ro = 2 ft~^Mpc, 7 = 1.8 and e = —1.2 proposed 
by Efstathiou (1995). Eor the spectral indices, we assume 
/3b ~ 1 — as the faint blue galaxies are predominantly late- 
type/lrr galaxies at increasing magnitude — and 13 x ~ 0.4 
appropriate for the spectrum of the unresolved XRB in both 
the soft and hard bands. Such a spectrum has been observed 
in absorbed AGN's and starburst galaxies and may therefore 
be appropriate for a large fraction of the faint blue galaxies. 
This parameter is actually barely significant when the oth- 
ers have been fixed. (Assuming /3x = 1 (bright AGN-like) 
changes the results below by less than 10%). The optical flux 
of a galaxy with magnitude m is Fb = lO""'^'™"'*' , where 
K = 16.6 for fluxes in units of lO^^h ergs s~^Mpc~^. The 
total optical emissivity due to galaxies (of all magnitudes) 
at redshift z verifies: 



^ These measurements have been carried out in the 0.4-2.4 keV 
band whereas / in defined in the 0.5-2 keV band. Using the spec- 
tral parameters quoted by the authors, we apply a mean spectral 
correction of Jo. 4— 2.4 = 1.3 /. 



Pb{z) 
4Trrf{z) 



dv{z) = dz 



10 



-0.4(r; 



' N{m, z)dm 



(32) 



Our adopted redshit distribution, arbitrarily extrapolated to 
B = 30 to account for 'all' the galaxies, yields Pb{z) ~ 3.8 x 
W^^{l + zy^^ h ergs s"^Mpc"^ out to z x 2. The correlation 
between optical and X-ray luminosities is approximately lin- 
ear for late-type galaxies, although very scattered (Eabbiano 
et al. 1988). Assuming Lx Lb at all redshifts and for all 
galaxies, then px(z) oc pb{z){1 + z)^^~^^ . In the present 
notation, this scenario corresponds to g = 1.3-|-/3_b — /3x = 2, 
and therefore V. = j — e + q — 5 = 0. The above values of K 
and ro imply (Eq. 27): 



Po 



2.7 X 10^** ih ergs s ^Mpc 



(33) 



for the volume emissivity of the galaxies at soft X-rays, and 
yield: 



log 



: log 



Po 



Pb{0) 



-3.1-Mog(z) 



(34) 



within the range of observed values for late-type galaxies 
(Eabbiano et al. 1988). This luminosity ratio implies that 
B > 18 galaxies have soft X-ray fluxes fainter than 10~^^ 
The source counts are shown in Eig. 3. 



ergs s ^cm^ 



EIGURE 3 

At the mean redshift of the sample z = 0.26, the comov- 
ing volume emissivity contributed by those X-ray sources 
which are clustered with the galaxies is: 



px{z) = po{l + zf ~ 4.2 X 10^** ih ergs s ^Mpc 



The integrated X-ray flux from 18 < B < 23 galaxies: 
Lx 

L^ " 



Ala 



(35) 



(36) 



dz{l + z 



10 



-0.4(r: 



N{m, z)dm 



represents ~ 22% of the XRB in the ROSAT band. The in- 
tegrated X-ray flux of 'all' (i.e. B < 30) galaxies depends 
on the maximum redshift of integration, and obviously on 
our arbitrary extrapolation of the redshift distribution above 
B ~ 24. The 'total' faint blue galaxy population, as mod- 
eled here, may contribute anywhere from 55 to 75 % of the 
soft unresolved XRB as Zmax increases from 2 to 5. Eigure 
4 shows the predicted correlation functions. The model is 
consistent with the upper limit of the auto-correlation func- 
tion of the background, although a slightly larger value for 
ro = 3 ft~^Mpc would improve the fit, as well as decrease 
the contribution of the catalogued galaxies to ~ 10%. 

EIGURE 4 



4.3 Further statistics 

Measuring the amplitude of the cross-correlation function 
between the diffuse XRB and increasingly faint galaxies can, 
in theory, allow us to significantly narrow down the param- 
eter space, and eventually derive fairly good estimates for 
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the volume emissivity and the clustering properties of the 
X-ray sources at increasing redshift. 

SAgt = H-,r'^^ {1 + zfrl-''N{m,z)dz Sm 

= ^(1 + z fri-^Nim, z)dz Sm (37) 

where /(K) is the function defined in Eq. 26. Figure 5a shows 
the normalized amplitude Swgb{m) = SAgb{m)/ IC{m) as a 
function of B— magnitude when K varies from —3 to 3. 

FIGURE 5a 

Measuring the cross-correlation function at zero-lag can, 
again in theory, yield the galaxy contribution to the back- 
ground as a function of magnitude {SIg{m)) in a straight- 
forward way. Equations 17, 18 and 24 imply: 

Srio{a) = a^5Ig{m) -\- a'''~'^ C~f5Agb{m) (38) 

where the cross-correlated cells are assumed to be square 
solid angles . Consequently we can hope to derive the X- 
ray contribution of a flux limited sample as well as the am- 
plitude of the angular correlation function, by measuring rjo 
for different cell sizes. In the magnitude range of the data: 

■no{a)=aAIg + a~''C.^A^. (39) 

Fig. 5b shows the lower and upper limits of the normalized 
zero-lag cross-correlation functions: wo = rjo I IN, when 
A/g/J varies from 10 to 100% of the soft XRB. The soUd 
line is our model prediction, i.e. l\lgjl ~ 22%. 

FIGURE 5b 



5 CONCLUSIONS 

This paper presents a formalism for the cross-correlation 
function between a diffuse background radiation, such as 
the XRB, and a population of discrete sources, such as the 
so-called faint blue galaxy population. By extending earlier 
modeling of the cross correlation of the XRB with galaxies 
at zero redshift, the present formalism allows to constrain 
the X-ray emissivity and clustering properties of these faint 
galaxies and other X-ray sources clustered with them. 
Cross correlation analyses of the hard (~ 2 — 10 keV) XRB 
with local {z < 0.03) optical and IRAS galaxies yielded 
po ~ 0.7-1.2 X lO^^ft ergs/Mpc^ (Lahav et al. 1993; Miyaji 
et al. 1994; Carrera et al. 1995). However, these results did 
not constrain the contribution to the XRB from high red- 
shift. Roche at al. (1995) have recently measured the cross- 
correlation function between faint galaxies in the magnitude 
range B = 18 — 23 and the soft (~ 1 keV) X-ray back- 
ground, but only provided a heuristic interpretation to their 
measurement. An application of our formalism to a specific 
model predicts that these 18 < B < 23 faint galaxies (cov- 
ering redshifts z < 0.5) make up ~ 22% of the 1 keV XRB. 
When extrapolated to fainter magnitudes, the faint blue 
galaxy population can account for most of the residual back- 
ground at soft energy. Our model satisfies all the available 



constraints on the 3 correlation functions (galaxy /galaxy, 
background/background and galaxy /background). However, 
we do not consider here the so-called 'spectral paradox', i.e. 
that the spectrum of the XRB (in both soft and hard bands) 
is not compatible with the spectrum of known sources. It re- 
mains an observational challenge (and a test for our model) 
to accurately measure the X-ray spectrum of the faint blue 
galaxies and other high redshift objects. 

While the idea that the faint blue galaxies could significantly 
contribute to the XRB seems plausible, we wish to empha- 
size several uncertainties in the current measurements and 
modeling, and to propose other strategies. The contribution 
A/g/J ~ 22% we have derived for the 18 < B < 23 galaxy 
population results from very specific assumptions. However 
appropriate these may be, different scenarios may lead to 
quite different conclusions. For example, a strong clustering 
model with ro = 5 ft~^Mpc (as observed for bright optical 
galaxies) and other previous assumptions unchanged, pre- 
dicts A/g// ~ 4%. As the contribution drops, the model re- 
mains consistent with the background auto-correlation func- 
tion. Inversely, very weak clustering with ro = 1 ft~^Mpc 
and e = 0, results in A/g// ~ 95% of the background. An- 
other option may be that the faint blue galaxies emit very 
negligible X-ray (A/g Si 0) but are strongly clustered with 
other sources of the XRB (e.g. very faint AGN's, or hot 
gas confined in poor galaxy groups), thus inducing a cross- 
correlation signal. 

On the observational side, it seems that important quanti- 
ties such as the mean soft XRB intensity and spectral in- 
dex, are still uncertain, in part due to the different ways 
'resolved' sources and the galactic component (strong below 
~ 1 keV) are removed in various analyses. In deriving the 
non zero-lag cross-correlation it is also crucial to carefully 
take into account the detector's point-spread-function, in 
order to make sure that self-correlation of objects does not 
contribute to the signal. Regarding data analysis, the zero- 
lag cross-correlation is in fact easier to measure, taking into 
account the beam smearing (cf. Miyaji et al. 1994), and more 
informative than the non zero-lag correlation. It would also 
be interesting to apply these cross-correlation techniques to 
higher redshift objects (e.g. quasars and clusters) and in 
harder X-ray energy bands where most of the XRB energy 
resides, and to relate it to the emissivity derived from the 
cross-correlation with local galaxies. As demonstrated in this 
paper, the cross-correlation technique is a powerful probe of 
the X-ray emissivity as a function of redshift. 
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where the average is over all cells. We assume the cells can 
be divided into elements SVk containing rik = or I galaxy 
(cf. Peebles 1980): 

N = J2"i' 

I = ^ ik, where ik = (pk/'inrl) SVk 
NI = ^ nkik + ^ rijik 



<NI> : 
< NX I> 

Therefore: 



4nrf 



■dV + 



nidV^i 



P2 



4nr? 



P2 



■dV2 



■[1 + i(ri2)]dVidV2 



Ip + Ic 



(A2) 



where rjp is the poisson term and rjc the clustering term 
defined as follows: 



dV 



P2 

4nr? 



■{(ri2)dyidy2. 



(A3) 



(A4) 



Pg is the background emissivity originating from the cata- 
logued galaxies themselves, and so rjp is the intensity con- 
tributed by these galaxies to the background (which may 
arise from just a fraction of them), rjc arises from the clus- 
tering of the galaxies with the background sources. The in- 
tegrals are over Vceii = ^ / dv(z). ri and r;^ are luminosity 
distances and ri2 is the physical separation between dV^i and 
dV2 . Distances and volumes are defined in Appendix B for 
different cosmologies. 



Al The Poisson term 

The background we are concerned with in this paper is the 
X-ray background. If the galaxies have been observed at 
frequency v and Lx{Li,) is the potential X-ray luminosity of 
a galaxy with luminosity Li,, the volume emissivity of the 
galaxies in X-ray may be written: 



APPENDIX A: ZERO-LAG 
CROSS-CORRELATION FUNCTION 

We start by deriving the zero-lag cross-correlation be- 
tween the galaxies and the background. The angular cross- 
correlation and auto-correlation functions follow quite sim- 
ply from the same principles. We divide the sky into cells of 
small solid angle w and calculate the following quantity: 



^a=< 5N5I >cells 

= < {N -N){I -T) > 
= < NI > - < N X I > 



(Al) 



Pg{z) = (I + Z) 



Lx{L,)<j){L,,z)dL,, (A5) 



where (f>(Li,,z) refers to the galaxy luminosity function at 
v and z, and jSx is their spectral index in the X-ray band. 
Lmtn{z) and Lmax{z) depend on the flux limits of the survey 
and the spectral index of the sources at i^. The mean X-ray 
intensity contributed by the catalogued galaxies is: 



AL 



4nrf(z) 



dv{z) 



and therefore: 



Ip 



(A6) 



(A7) 
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A2 The clustering term 

If we assume that sources at significantly different redshifts 
are uncorrelated ( { = when |r2 — »"i | >> 0), we may do 
the following change of variables, in comoving coordinates: 



Xl={l+ z)Xc 

dyidy2 = F^X^dXcdQidQ2du 



(A8) 



where F ■ 



dV/di, 



(=1 in flat space), as defined in Appendix 



B. The clustering term becomes: 
1c= [ f n(z)^^i(rr2,z)dVrdV2 



(A9) 



"i^)^^^cF^dXc I I dQidQ2 / iiri2,z)du 



In the small angle approximation, the physical separation 
between the 2 sources is: 



A2.1 power spectrum formalism 

^ 1 /"°° 7^/7 ^ sin(A;r/a) , 2 



(AlO) 



(All) 



where k is the comoving wavenumber and a = [1 -\- z) ^ is 
the scale factor. The integral over u is: 

-t- CO {'OO 

i{ri2,z) du = 2 I i{ri2,z)du 



A2.2 power-law approximation 
{(r,^) = (l + ^)-(^+^' {- 



(A16) 



The integral over u becomes: 

-t- CO 

^(ri2,z) du 



-\- CO 

(l+^)-(3 + ^)+7,7 / [u^F^ +Xle^)-''l^du 



{l + z)-'-^+^'>+-W-iF-'H^{XJ) 



1-7 



2^,-7/2 _ ^ 12^^ I 2 
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r(i) 



(A17) 



(A18) 



For a square beam with size lo = a^ , the integral over d can 
be written (Totsuji & Kihara 1969): 



dQidQ26i^"'' = C-,a^"'' 



Finally, 

ric = a^~'^ H -.fC-.fr g x 



(A19) 



(A20) 



"'^'^l^Xflzj'^^ + ^)-(^+^'+-'Xr-'(^)F(^)dXe(^) 



For 7 = 1.8, if-, = 3.68 and C\ = 2.25. 



CO /"CO 



,sin(A;r/a) ,2 , 
P(fc,g) k dk du 

{kr/a) 



(A12) 



The inside bit is: 

sin(A;r/a) 1"^ sm 

(fcr/a) " " (p„2f2 _^px2e2)i/2 



sm 

_^PX26l2)l/2 



F-i|jo(A;Xee) 



.d„=/ ^d« 



(A13) 



where Jo is a Bessel function defined as: 

, , , 2 f°° sin(t^+xy^ , 2 P . , , , , 
Jo(x) = — I — j-^j^ — ; — /o — dt = — / sm 2;(cosh u)du. 



(A14) 



■ (t2+2;2)l/2 

Inserting this into rjc, we can write: 

/• CO 

ric = P{k, z)kgo{k)dk 
Jo 

where the kernel go is defined as: 

~g,(k) = 1- jn(z)^^^Xt(z)F(z)dX.(z) 

X dnidn2-Jo{kXc0) 



(A15) 



Similar calculations have been presented by Baugh & Ffs- 
tathiou (1994) for the galaxy auto-correlation. 



A3 Angular correlation functions 

The background/background and galaxy/galaxy auto- 
correlation functions and the galaxy /background cross- 
correlation function at angular separation 9 are respectively 
defined as: 

r]bb{0) =< Sle'SIe' + e > = < h'h' + e > -P 

rigg(ff) =< SNe.SNe. + e > = < Ne-Ne.+e > -N^ (A21) 

rigb(d) =< SNe'SIe' + e > = < Ne'h' + e > -NI 

There is no Poisson terms this time since the cells do not 
overlap (^ > 0), simply: 

nM0) = lJ;^J^mliiri2,z)dV2 

1gg(S) = J^nin\dViy i(ri2,z)dV2 (A22) 
rj,b{0) = Jlr^{z)^£|i§L.dV^Jarl2,z)dV2 

With the same change of variables as in the previous section 
(Fq. A8), we can write: 

r^tbie) = C P{k,z)kgbb{ke)dk 

= Jo°° Pik, z)kg,,{k0)dk (A23) 
ngb{0) = J^°° P{k,z)kggb{k0)dk 

where 



gbb(k0) = ^l[f^) Jo(kXJ)Xt 
ggg{k0) = i^l n'{z)Jo{kX,e)XtFdX, (A24) 
ggbikff) = 2. j^n{z)^j,(kXJ)XtFdX, 



or else, in the power-law approximation: 
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rj,,{0) = H,r-<0'-^ I (i^) (1 + zyxt^FdX, 
^gg(e) = H^r-^e'--' n'(z)(l + zfXt^FdX, (A25) 
rj,,{0) = H,r-<0'-^ j\{z)£^{l + zyxt^FdX, 

where p = 7 — (3 + e) and again Xc{z) is the comoving 
distance coordinate rdz) and Xi{z) is the corresponding lu- 
minosity distance defined in Appendix B. Here we have as- 
sumed the size of the cells to be infinitesimal dco, so that the 
previous integrals dflidfl2Jo{kXcO) or, in the power- 

law approximation, J J dQidQ2^^~"', over the cell size are 
replaced in the above equations by dLOidL02Jo{kXcO) and 
dci)idci)2^^~"' respectively. All the above correlation functions 
are thus implicitly defined per unit square solid angle. 



APPENDIX B: COSMOLOGY 

In the following, fc is the comoving radial coordinate in units 
of Mpc and dv = dV/dco the volume element per unit solid 
angle in units of Mpc^sr"^ . F is a curvature function defined 
as F = dv/r^dfc (=1 in flat space). The luminosity distance 
in Mpc is defined as ri{z) = (1 -|- z)rc{z). 



Bl Qo = 1 and A = : Flat (Einstein-de Sitter) 

r4z) = {2c/Ho)[l-il + z)-'''] (A26) 
(c/Ho) 



dz rldz ~ (1 + ^)3/2 



(A27) 



B2 Qo < 1 and A = : Open 

(Friedmann— Robertson— Walker) 

= (2o/H., "°'-^'"';|;K;^y°"'"-" ,A..8) 



B3 Qo < 1 and A = 3(1 - 9.o)Hl : Flat with a 

Cosmological Constant (Friedmann— LemaTtre) 

rJz) = (c/Ho) / ^ (A30) 

F(z)^(z) = 4^ = — (A31) 
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